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ABSTRACT 

We present observations of high temperature emission in the core of a solar active 
region using instruments on Hinode and SDO. These multi- instrument observations al- 
low us to determine the distribution of plasma temperatures and follow the evolution 
of emission at different temperatures. We find that at the apex of the high tempera- 
ture loops the emission measure distribution is strongly peaked near 4 MK and falls off 
sharply at both higher and lower temperatures. Perhaps most significantly, the emission 
measure at 0.5 MK is reduced by more than two orders of magnitude from the peak at 
4MK. We also find that the temporal evolution in broad-band soft X-ray images is rel- 
atively constant over about 6 hours of observing. Observations in the cooler SDO / AIA 
bandpasses generally do not show cooling loops in the core of the active region, con- 
sistent with the steady emission observed at high temperatures. These observations 
suggest that the high temperature loops observed in the core of an active region are 
close to equilibrium. We find that it is possible to reproduce the relative intensities of 
high temperature emission lines with a simple, high-frequency heating scenario where 
heating events occur on time scales much less than a cooling time. In contrast, low- 
frequency heating scenarios, which are commonly invoked to describe nanoflare models 
of coronal heating, do not reproduce the relative intensities of high temperature emission 
lines and predict low-temperature emission that is approximately an order of magnitude 
too large. We also present an initial look at images from the SDO / AIA 94 A channel, 
which is sensitive to Fe XVIII. 

Subject headings: Sun: corona 



1. Introduction 

Understanding the origin of high temperature plasma in the solar corona is one of the central 
problems in solar physics. The nanoflare concept represents one of the more popular theories for de- 



scribing how energy stored in the Sun's magnetic field is converted into thermal energy (e.g. JParker 
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19721 . Il983l ) . In nanoflare models turbulent photospheric motions drive the twisting and braiding of 
the magnetic field, which leads to the release of energy on small spatial scales as the accumulated 
topological complexity is dissipated by magnetic reconnection. Nanoflares are often modeled as 
impulsive heating events where the plasma is heated to very high temperatures (~ 10 MK) and 
then cools down through a combination of conduction, enthalpy flux, and radiation without be- 
ing reheated. Furthermore, since these heating events are likely to occur on very small spatial 



scales an observed loop is assumed to be composed of ma ny unresolved strands (e.g., ICargill 



1994 



Klimchuk fc Cargilll l200ll : iPatsourakos &: Klimchu k 2006). This framework implies that nanoflare 
heated loops should have co-spatial hot and cool emission. 

Observations at relatively cool coro nal t emper a tures , however, have have cast doubt on this 
heating scenario. lAntiochos et~aD (|2003h and lNittal (hoooh . for example, have argued that cooling 
loops are often not observed in the core of an active region. This position has been supported by 
recent observations with the EUV Imag i ng Spectrograph (EIS ) and the X-ray Telescope (XRT) on 
the Hinode mission ( Kosugi et al.1 120071 ). I Warren et al.1 ( 20101 ) found no evidence for emission near 
1 MK that was spatially correlated with emission at higher temperatures (3-5 MK) in the core of 
an active region. In this region the high temperat ure emission measured wi th XRT was observed to 
be relatively constant over many hours. Similarly. iBrooks Warrenl (|2Q09l ) have found no evidence 
for strong variability in the intensities and Doppler shifts measured in the moss, the footpoints of 
the high temperature loops. These results suggest that high temperature loops are generally heated 
on very short time scales, much shorter than a characteristic cooling time, and do not cool down 
to lower temperatures. 

Previous studies on the properties of high temperature emission in the core of an active region 
have failed to adequately address a fundamental question: what is the distribution of temperatures 
in the core of an active region? The EIS spectral range contains emission lines from Ca XlV-Ca XVII 



which provide exce llent coverage of the critical 3-5 MK temperature range (see iDel Zanna 



2008 



Warren et al.ll2008l ). Observations of these lines in combination with other emission lines observed 
with EIS and observations with the thick XRT filters allow for the emission measure distribution to 
be computed over a very wide range of temperatures. The properties of the loop apexes along the 
line of sight are of particular importance. The dispersion in the temperature constrains the rate at 
which energy is input into high temperature loops and such information is critical for motivating 
physical models of this emission. 

In this paper we present the analysis of new observations in the core of an active region 
taken after the launch of the Solar Dynamics Observatory (SDO). These observations feature 
spectroscopic observations from EIS and soft X-ray images from XRT on Hinode. These data 
show that the emission measure distribution for the loop apex is peaked at temperatures near 
4MK and falls off sharply at both higher and lower temperatures. At temperatures near 0.5 MK 
the emission measure is more than two orders of magnitude below the value at 4MK. Simple 
hydrodynamic modeling suggests that these properties are inconsistent with low frequency nanoflare 
heating models. Such models predict low temperature emission that is approximately an order of 
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Fig. 1. — AIA and XRT observations of active region 11089 taken on 2010 July 23 near 15 UT. 
The field of view is 384" x 384", which is motivated by the size of the XRT images. The dotted line 
indicates the extent of the EIS raster. The small box indicates the region between the moss used 
to calculate the emission measure with EIS and XRT, and represents the physical conditions at 
the loop apexes along the line of sight. This region was chosen because it minimizes contamination 
from the moss. The electronic version of the manuscript includes movies of the AIA 171, 335, 94, 
and 131 A channels. 

magnitude larger than what is observed. 

The high cadence (~ 12s), high spatial resolution (0.6" pixels), narrow band images provided 
by the Atmospheric Imaging Assembly (AIA) on SDO provide additional information on the evolu- 
tion of plasma at a wide range of temperatures. Images from channels dominated by 1 MK emission 
show only a few, isolated cooling loops in the moss regions, consistent with previous results. In 
this paper we also present some initial observations from the AIA 94 A channel, which is sensitive 
to the Fe XVIII 93.94 A line formed at about 7MK. In small flares and microflares the expected 
progression from high temperature to low temperature emission is observed. In more quiescent 
areas of the active region core, however, the emission appears to be relatively constant, potentially 
consistent with what is observed in XRT. 
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Fig. 2. — EIS observations of AR 11089 in various emission lines. The field of view is 120" x 384". 
The small box indicates the region between the moss used to calculate the active region "core" 
emission measure. The size of this region is 10" x 15" or 75 spatial pixels. 
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Fig. 3. — XRT observations of AR 11089. On the left is a frame from the observing sequence. On 
the right are light curves from three different locations in the core of the active region. Intensities 
are in DN s _1 per pixel. The dashed horizontal line is the median intensity for the period shown. 
The dotted horizontal lines are ±20% of the median. The solid vertical line corresponds to the time 
of the image. The emission in the core of the active region is relatively constant, with a variability 
of 20% or less. Transient events with relatively short lifetimes (< 900 s) are also observed. The 
electronic version of the manuscript includes a movie of these data. 



2. Observations and Analysis 

In this paper we consider EIS, XRT, and AIA observations of NOAA active region 11089 taken 
on 2010 July 23. This section gives a brief overview of the various instruments, the available data, 
and the calculation of the differential emission measure distribution. 

The AIA on SDO consists of 4 independent normal incidence, multi-layer telescopes, each 
with 2 channels. Each channel images the full Sun. The AIA coronal channels and their dominant 
components are as follows: 94 A (Fe X and Fe XVIII), 131 A (Fe VIII and Fe XX, Fe XXI, and 
Fe XXIII), 171 A (Fe IX and Fe X), 193 A (Fe XII), 211 A (Fe XIV), 335 A (Fe XVI). A detaile d 
discussion of the predicted contributions to each AIA channel is presented by lO'Dwyer et al.l (|2010l ). 
The relatively cool components of t he 94 and 131 A channels have been iden tified in high spectral 



resolution stellar observations (e.g.. iRaassen et al. 



2002 



Liang Sz Zhaoll2010l ). There is also a He II 



304 A channel and a UV channel. Figure Q] shows sample AIA images from this active region 
region taken near 15 UT. For this work we use the level 1.5 cutout data that has been re-scaled 
and co-aligned to a common coordinate system. We note that to precisely co-align the data in 
the various channels we cross-correlated a subset of 20 full-disk images taken from the period of 
interest. These offsets were generally on the order of 1-2 pixels in each direction. De-rotation of 
the data was achieved by using established solar rotation rates. As is evident in the movies, there 
is very little jitter in the data. 
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The XRT on Hinode (jGolub et al.l 120071 ) is a grazing incidence, soft X-ray telescope. Temper- 
ature discrimination is achieved through the use of focal plane filters. Because XRT can observe 
the Sun at short wavelengths, XRT can observe high temperature solar plasma very efficiently. 
For the observations considered here, XRT took images in the thin Ti-poly filter at a cadence of 
about 120 s, with Al-thick context images taken approximately every 90 minutes. The observing 
was periodically interrupted by seasonal eclipses that occur annually during the May - August 
period. Representative XRT Ti-poly and Al-thick images are shown in Figured! The XRT images 
have been co-aligned with respect to the first image in the sequence. The XRT and AIA data were 
co-aligned by cross-correlating the Ti-poly and 335 A images. This analysis incorporates the latest 
available calibration for XRT (jNarukage et al.ll2010l ). 



The EIS instrument on Hinode ( Culhane et al. 



2007 



Korendyke et al.ll2006l ) is a high spatial 



and spectral resolution imaging spectrograph. EIS observes two wavelength ranges, 171-212 A and 
245-291 A, with a spectral resolution of about 22 mA and a spatial resolution of about 1" per pixel. 
Solar images can be made by stepping the slit over a region of the Sun. Telemetry constraints 
generally limit the spatial and spectral coverage of an observation. For the observations considered 
here the 1" slit was stepped over the core of the active region using 60 2" steps and 512" of the slit 
height was read out. The field of view for the EIS raster that was run between 14:32 and 15:34 UT 
on 2010 July 23 is indicated in Figure [TJ 

The use of wide 2" steps in the EIS raster degrades the spatial resolution somewhat, but allows 
for the raster to be completed in between Hinode eclipses even while using relatively long (60 s) 
exposure time. The use of wide steps also reduces telemetry usage and allows for a very extensive 
line list to be telemetered to the ground. From these data we have determined intensities for 39 
spectral lines. Images for some of these emission lines are shown in Figure [2j Almost all of the line 
intensities can be determined by simple Gaussian fits. The Ca XVII 192.858 A line is complicated by 
blending with Fe XI 192.813 A and a complex of O V lines, inc luding O V 192. 906 A. The intensities 
for these lines are determined using the method outlined by iKo et al.l (|2009) , where the intensity 
of the Fe XI 192.813 A line is inferred from Fe XI 188.219 A and the remaining part of the profile 
is fit with multiple Gaussians. The resulting Ca XVII raster is consistent with Fe XVII 254.87 A, 
indicating that the deconvolution is approximately correct. 

Collectively, these data suggest a relatively simple morphology to the active region. The 
EIS emission lines formed at high temperatures, Fe XVI (2.82 MK), Ca XIV (3.55 MK), Ca XV 
(4.47 MK), CaXVI (5.01 MK), and Ca XVII (5.62 MK), and the XRT Ti-poly and Al-thick images 
show relatively short, hot loops that connect across the neutral line running through the middle 
of the active region. The temperatur es given here correspond to the peaks in the CHIANTI 6.0 
ionization fractions (jDere et al.ll2009l ). Below these hot loops are regions of intense "moss," the 
footp qints of the hot l o ops that are evident in emission lines formed a t temperatures near 1 MK 



(e.g., IPeres et al J 1 1994 : iBerger et al.lll999l ; Fletcher fe de Pontieu 



mentioned earlier, the relatively cool fan-like loops (e.g.. ISchriiver et al 



1999 



Martens et al.1 120001 ). As 



1999) appear to be largely 



absent from the core of the active region. The emission at these temperatures, however, is not 
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zero. There is certainly diffuse million degree coronal emission lying above the active region. The 
AIA movies also indicate that there is some million degree emission associated with the very cool, 
and highly dynamic filamentdike material in the core of the active region that connects across the 
neutral line at very low heights. 

As mentioned previously, the primary strength of the XRT is the ability to observe high 
temperature emission over a wide field of view at relatively high cadence. This is illustrated in 
Figure El where light curves from several positions in the core of the active region are shown. Each 
light curve is from a sin gle pixel, there is no spatial averaging. Consistent with earlier analysis 
(e.g.. IWarren et al.ll2010l ). these light curves show relatively constant emission. Typical intensities 



generally lie within a range of ±20% over many hours. Transient brightenings are observed, but 
they are relatively short lived, with lifetimes of 900 s or less. 

The AIA movies, which are included in the electronic version of this manuscript, show that the 
moss is consistent with the relatively constant emission observed in XRT. The AIA 131 and 171 A 
movies show clear evidence for many fine coronal loops in the active region. Few of these loops, 
however, connect to elements of the moss. There is evidence for intensity fluctua tions in the moss 



itself , but it is likely that this is related to spicular activity that obscures the moss (jde Pontieu et al 



1999). 



The primary objective of this study is to determine the distribution of temperatures in the 
hot loops in the core of the active region. We are particularly interested in the conditions at 
the loop apexes along the line of sight that avoids contamination from the moss emission. Such 
measurements are important for determining the physical conditions in the high temperature loops. 
To accomplish this we have determined spatially averaged EIS line profiles in the small "inter-moss" 
region indicated in Figures Q] and [2j The size of this region is 10" x 15" or 75 spatial pixels when 
we account for the 2" steps. This region has relatively little contamination from the moss or from 
overlying loops. The spatial averaging diminishes the uncertainties for the intensities of some of 
the weaker lines and allows for a more complete analysis. The intensities for these line profiles are 
computed using the same assumptions as the spatially resolved intensities shown in the rasters. 
The intensities from this region are given in Tab led! The selection of lines is motivated by p revious 



emission measure calculations in the quiet Sun ( Brooks et al 



2009; 



Warren Sz Brooks 2009) and in 



m 



flares (jDel Zannall2008l ; IWarren et al.ll2008l ). This previous work has identified emission lines that 
are mutually con sistent. Note that F e XVII 254.87 A has been left out because of uncertainties in 
the atomic data (jWarren et al.l l2008). 

To compute the differen tial emission measure we use t he Monte Carlo Markov Chain (MCMC) 
emission measure algorithm (IKashyap &; Drakdll998l . l2000l ) distributed with the PINTof ALE spectral 
analysis package. This algorithm has the advantage of not assuming a functional form for the 
differential emission measure. The MCMC algorithm also provides for estimates of the error in 
the DEM by calculating the emission measure using perturbed values for the intensities. The 
algorithm assumes the uncertainties in the intensities are uncorrelated so that systematic errors in 
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Log Temperature (K) 

Fig. 4. — The emission measure distribution for the core of the active region. The DEM was com- 
puted using the "inter-moss" intensities from EIS and XRT and represents the physical conditions 
at the loop-tops along the line of sight. The solid red line is the best-fit DEM and the dotted 
black lines are from Monte Carlo calculations using perturbed values of the observed intensities. 
250 Monte Carlo simulations of the emission measure were performed. This gives a statistically 
plausible range for the emission measure in each temperature bin. The solid lines are the emission 
measure loci curves described in the text and color-coded according to element. The emission mea- 
sure distribution is strongly peaked at logT = 6.6 and falls off sharply for both higher and lower 
temperatures. 

the calibration, which could depend on the wavelength, or in the atomic data, which could vary by 
ion, are not accounted for. 

In its current implementation the MCMC algorithm does not allow for the density to be a 
free parameter. In principal the Ca XV 181.900/200.972 A ratio should provide information on the 
density at temperatures near 4MK. The weaker component of this line pair, 181.900 A, however, is 
not observed in these data. At the densities we expect (< 10 10 cm' 3 ) the intensity of the 181.900 A 
line is less than 0.1 times that of the 200.972 A. This implies less than 20 DN would be detected in 
this 60 s observation. Monte Carlo simulations indicate that approximately 100 counts are needed 
to compute a line intensity. In the absence of alternatives we use the Fe XIII 202.836/202.044 A 
ratio, which yields a density of logn e = 9.45. We will discuss this result later in the paper. 
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The relationship between the emissivities, the differential emission measure, and the observed 
intensities is the usual expression 

h = -^ J ex(n e ,T e )Z(T e )dT e , (1) 
where e\(n e ,T e ) is t he emissivity compu ted with the CHIANTI 6.0.1 atomic databas e assuming 



coron al abundances (IFeldman et al.lll992l ) and the new CHIANTI ionization fractions (jDere et al 



2009). The function f(T e ) = n 2 e ds/dT e is the differential emission measure distribution inferred 



from the intensities. Here we also show the emission measure loci computed from 

£,loci(T e ) = / tTT> (2) 

e\(n e ,T e ) 

which indicates the temperature range where the various lines are sensitive. Note that to aid in 
the comparisons with the em loci we will always plot the DEM multiplied by the temperature bin, 

Z(T e )dT e , (3) 

and we refer to this as the emission measure distribution. 

The DEM computed from the MCMC method is shown in Figure 01 The intensities computed 
from the DEM, which are also given in Table [H are generally consistent with the observations to 
within 25%, although there are some significant discrepancies. We note that the Mg lines, which 
lie at the low end of the temperature range, are not fully consistent with Si VII. We conjecture that 
this is due to uncertainties in the abundances. Since the intensity for Si X is consistent with the 
other Fe lines formed at a similar temperature, we assume that this discrepancy is due to the Mg 
abundance. To compute the emission measure we multiply the coronal Mg abundance by a factor 
of 1.7 to bring Mg VII 280.737 A into agreement with Si VII 275.368 A. In this analysis we have also 
included two lines from S, which is a low first ionization potential element. The consistency between 
the S and Fe lines in this DEM calculation suggests that the assumption of coronal abundances 
is generally valid, at least in a relative sense. Among the Fe lines the most significant problem is 
for Fe XIII. This is perplexing since there are no such problems evident in the analysis of quiet 
Sun data. This issue is unresolved here, but the emission measure is meant to be only a crude 
representation of the temperature distribution in these core active region loops. We will consider 
comparisons between some simple heating models and the observed intensities in the next section. 

The most significant result from this calculation is the fact that the emission measure is strongly 
peaked at about logT e = 6.6 (4MK) and falls off sharply at both higher and lower temperatures. 
The peak emission measure is approximately 10 28 cm -5 while the emission measure at logT e = 5.8, 
the temperature of formation for Si VII, is about a factor of 400 less. Because of the broad 
nature of the Ca XVII and XRT Al-thick responses, the decline in the emission measure at higher 
temperatures is less certain. The Monte Carlo simulations suggest a fairly broad range of possible 
emission measures at the higher temperatures, but the vast majority of them lie at the lower end 
of the range. 
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Numerous a c tive region emission measure distribution s have been publish e d in the literature 



(e.g., Dere 



1982 



201ll : lTripathi et al 



2011 



Dere Sz Mason! Il993l : iBrosius et al.l Il996l ; IWarren et al.l 1200 ll ; IWinebarger et al 



Our emission measure distribution has significantly steeper slopes away 
from the maximum than most previous results. We attribute this to the fact that most previous 
work has used intensities derived from averages over large field of view. Our work, i n contrast, has 



focuse d on a very small region tha t represents the conditions near the loop apex. iTripathi et al 



(|201ll ) and IWinebarger et al.l (|2011l ) have also considered a small "inter- moss" region and find both 
shallow and steep slopes from the peak emission measure to lower temperatures. It is clear that 
a systematic study is required to understand the range of possible emission measure distributions 
near the loop apex. 

In computing the intensities in the active region core or inter-moss region we have not per- 
formed any background subtraction. The intensities of the high temperature lines outside of the 
core are negligible. For the cooler emission lines it is not clear which region to take as the back- 
ground. The intensities for these lines in the dimmest areas of the active region are typically 10% 
of the core intensity and close to the core the intensities are comparable. This suggests a very 
large range of possibilities. In all cases the background subtracted intensities would be lower than 
the values that we have used. This would lead to even steeper declines in the emission measure 
distributions at low temperatures. At low temperatures, the emission measure shown in Figure U] 
is an upper bound. 



3. Modeling 



Early work on the modeling of high te mperature loops observed at soft X-ray waveleng ths of- 
ten focused on steady heating models (e.g-. lKano k, Tsunetalll995l : lPorter Klimchuk!ll995l). Such 



analy sis has also yielded volumetric filling factors that were generally less than 1 (IPorter &: Klimchuk 



19951 ). indicating that the observed loops were unresolved. Such observations could not exclude the 
possibility that the observed emission only appeared steady, but was actually made up of many 
unresolved loops that are evolving. The broad temperature coverage of the combined EIS and XRT 
observations allows us to make much more detailed comparisons with the predictions of different 
heating scenarios. 

In this section we will make quantitative comparisons between the observations and solutions 
to the hydrodynamic loop equations, which describe the evolution of mass, momentum, and energy 
in a flux tube. We will consider low-frequency heating, where the time between heating events is 
long relative to the cooling time, and high-frequency heating, where the loop is reheated before it 
has time to cool. The goal is to consider some very simple cases that will lay the foundation for more 



(e-g- 


Schriiver et al. 


2004; 1 


Warren & Winebarger 


2007 




Winebarger et al. 


2008; 


2008; 


Winebarger et al. 


2011 


.) is beyond the scope of this work. 
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Fig. 5. — An example of high-frequency heating that keeps a loop close to equilibrium. Shown 
here are the apex temperature and density (top panels) in the loop and the apex intensities for 
four high temperature emission lines (bottom panels). For the line emission the average intensity 
is also indicated. The initial atmosphere for this simulation is very tenuous (0.4 MK initial apex 
temperature) and the density converges to the asymptotic value slowly so the averaging is done 
over the final 8,000 s of the simulation. 

The loop length is an important parameter in any hydrodynamic simulation. The AIA images 
suggest a distance of approximately 50 Mm between the middle of the moss regions in this active 
region. This implies a total loop length (L) of approximately 75 Mm for a semi-circular loop that 
is perpendicular to the solar surface. For simplicity we will consider a single loop length in the 
simulations. A radiative loss rate as a function of temperature must also be specified. The DEM 
calculation presented in the previous section showed that the observed intensities for S X and S XIII 
are consistent with coronal abundances so we will assume a radiative loss curve computed using 



-12- 



the coronal abundances of lFeldman et al.l ()1992l ) 



To solve the time-dependent hydro dynamic loop equations we use the NRL Solar Flux Tube 
Model (SOLFTM). We adopt many of the same parameters and assumptions that were used in 
previous simulations with thi s code and we refer the reader to the earli er papers for additional details 
on the numerical model (e.g., iMariska et al.lll989l : IWarren et al.ll2003l ). To establish a characteristic 
time scale for the cooling of a high temperature loop with a total length of 75 Mm and a coronal 
composition, we begin by considering an initial simulation of a loop in equilibrium with an apex 
temperature of 4MK, the peak temperature in the DEM. If the loop is allowed to cool without any 
additional heating the apex temperature reaches 0.5 MK in approximately 1200 s, we will use r c to 
represent this cooling time. 

To model emission that persists we need to specify the volumetric heating rate as a function 
of time. We consider a series of step-function heating events that have a magnitude e, a duration 
5, and an occurrence rate r. We will assume that 5 < r for these simulations. We will also assume 
that the heating is uniform over the loop length. The mean heating rate is simply 

-- til 

T 

For a given loop length and radiative loss curve the parameter e determines, at least approxi- 
mately, the average apex temperature for the loop. For this work we will assume e = 8.31 x 
10~ 3 erg cm -3 s _1 , which is the heating rate required to keep the loop in equilibrium with an apex 
temperature of approximately 4MK, and we will vary e, 5, and r subject to this constraint. 



(4) 



An estimate of the event occurrence rate was given by ICargiill (119941 ) using 

Q 



T ~ 



AA h : 



(5) 



where Q is the average event energy, A is the radiative loss rate of the corona, and Ah is the spatial 
scale for coronal loops. For Q ~ 10 24 erg, ~ 10 cm 2 , and A ~ 10 7 erg cm -2 s -1 the time scale 
is r ~ 1000 s. Since these parameters are simply estimates there is a considerable range of possible 
values for r. With the application of the appropriate geometrical factors Equation [J] can be recast 
in a form similar to Equation [5) Our preference is to use the volumetric heating rate since it is the 
input to the hydrodynamic loop equations. 

For high-frequency heating (r <C r c ) we expect the loop to be close to equilibrium. To illustrate 
an example of such a loop we consider the parameters r = 150 s, 5 = 13 s, and e = 9.28 x 
10~ 2 erg cm -3 s _1 . To simplify the discussion we average the solution over the top 50% of the loop 
length at each time step to determine the apex density and temperature as a function of time. We 
then compute the time-dependent intensity of each emission line of interest using 



h 



-^e x (n e ,T e )n 2 e ds. 



(6) 



We infer the path length (ds) from the observed Ca XV 200.972 A intensity and use this value in 
calculating the intensities for all of the other lines. For this example the path length is 12.3 Mm. 
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Finally, to approximate the superposition of many sub-resolutio n strands in various stages of 



heating and cooling we com pute the time-averaged intensity (e.g., iPatsourakos &: Klimchuk 



2006 



Warren &: Winebargerll2007l ). The resulting apex temperatures and densities, are shown in Figure [5] 
along with the intensities for the high temperature Ca lines. All of the simulated and observed 
intensities are given in Tabled 

For this set of parameters the high temperature emission lines (above log T = 6.45) are gener- 
ally well produced by the model. Since heating events occur very frequently the the temperatures 
and densities in the loop are never too far from their average values (log T = 6.60 and log n e = 9.64), 
and the loop never cools. These simulation results are from a family of simulations where the heat- 
ing duration 5 and heating rate e were varied. For larger values of delta the departures from the 
average temperature is smaller and the intensity in Ca XVII becomes much smaller than what 
is observed. Similarly, for smaller values of delta the intensities in Ca XVII are larger and also 
inconsistent with what is observed. 

The intensities for the lower temperature emission are not reproduced by this model. Since 
the loop never cools, the modeled intensities at low temperatures are close to zero. One possible 
interpretation is that the relatively weak cool emission is unrelated to the high temperature active 
region loops and comes from either very low-lying loops or very long loops that extend over the 
active region. 

For r > r c the loop will evolve over a large range of temperatures and will be far from 
equilibrium. To illustrate an example of the low-frequency heating that is usually assumed in nano- 
flare modeling we consider the parameters r = 1800s, 5 = 67s, and e = 2.22 x 10 _1 erg cm -1 s _1 . 
The apex temperatures and densities from this simulation are shown in Figure [6) These parameters 
give approximately the same mean temperature and density (logT = 6.54 and logn e = 9.62) as in 
the previous case. In addition to the high temperature Ca lines we also show the evolution of Fe XII 
195.119 A and Si VII 275.368 A in Fi gure [H The path length inferred from the Ca XV intensity is 
25.2 Mm. All of the simulated and observed intensities are given in Table El 

In this case the long time between heating events allows the loop to cool and there is significant 
emission for lines formed at low temperatures. These intensities, however, are much larger than 
what is observed. For Si VII 275.368 A, for example, the modeled intensity is almost a factor of 
10 times too large. The difficulty for the low-frequency heating scenario is the s low draining of 



the d ensity from the loop. Previous wor k has suggested an n r ~ y/ T^. relationship (jJakimiec et al 



19921 : see also the more recent work by iBradshaw fc Cargilll 120101 which refines this scaling law) . 



This leads to relatively large densities and thus large intensities at low temperatures. 

The low-frequency heating case also has difficulties with the high temperature emission. The 
modeled intensities for Ca XVII 192.858 A and XRT Al-thick are larger than what is observed. 
This simulation is from a family of simulations with different values for 5 and e. For the cases we 
considered (5 < 100 s) the intensities for this emission was always about a factor of 3 greater than 
what is observed. This seems to be an inevitable consequence of low frequency heating. To produce 
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relatively high densities with short bursts of heating the magnitude must be large, which leads to 
high temperatures and large Ca XVII 192.858 A and XRT Al-thick intensities. 

As mentioned previously, the power law index that we measure in the emission measure dis- 
tribution between the peak near logT = 6.6 and lower temperatures is steeper than most previous 
measurements and steeper than is ex pected from s imple heating models, which can have a emission 



measure that scales like ~ T 3 / 2 (e.g.. |jordanlll976l ). This is likely a result of averaging over a small 



region near the loop apexes in the observations. To demonstrate this we have calculated the emis- 
sion measure distributions, essentially histograms of n^ds, from both the high- and low-frequency 
heating simulations. These distributions, which are shown in Figure [71 clearly illustrate the effect 
of averaging. The power-law index (6 in EM ~ T b ) for the entire loop is smaller (6 < 1.5) than for 
the loop apex (b > 2). 

One objection that could be raised regarding the emission measure distribution is that the 
box used to compute the average intensities is small and these values are not representative. The 
contrast between the moss and inter-moss regions at low temperatures, however, is small. For 
Si VII 275.368 A, for example, the intensities only increase by a factor of about 2 in the moss. It is 
clear that the large discrepancy between the intensities computed from the low-frequency heating 
scenario would not be changed significantly by modifying the region used for the averaging. 

Since the rasters for the lowest temperature lines are relatively noisy, it is tempting to think 
that these lines are simply too faint to measure reliably. As we illustrate in Figure El however, 
the emissivity for Si VII 275.368 A, which is one of the primary EIS lines at low temperatures is 
actually higher than the emissivity for the high temperature Fe XVI 262.984 A line, which is easily 
observed in the core of the active region. The effective area for EIS at these wavelengths is also 
very similar. Absent some fundamental flaw in the hydrodynamic models, it appears that any hot 
plasma evident in the Fe XVI 262.984 A line would be easily observed in Si VII 275.368 A as it 
cooled. 

The emissivities for the hot lines shown in Figure [8] suggest similar inconsistencies between the 
observations and low- frequency, impulsive heating scenarios. The emissivity for Ca XVII 192.858 A 
is higher than that for Ca XIV 193.874 A, suggesting that loops that have been heated up to very 
high temperatures (~ 10 MK) and are cooling would show comparable intensities in these lines. 
Instead, the data show that the Ca XVII emission is weaker than the emission from Ca XIV. Again, 
the observations are difficult to reconcile with traditional hydrodynamic models of impulsively 
heated loops. 



4. AIA 

As we have already seen in Figure[H AIA's combination of high cadence, high spatial resolution, 
and broad temperature coverage provide new ways to explore the temperature structure of an active 
region. Of particular interest to studies of high temperature active region loops is the 94 A channel, 
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which contains the Fe XVIII 93.94 A formed at about 7MK. Also of interest is the 131 A channel, 
which contains Fe XX, Fe XXI, and Fe XXII lines. Another channel not previously flown is the 335 A 
channel, which images Fe XVI 335.4 A. In Fi gure [8] we show the isothermal temperature response 
curves for the AIA 171, 335, 94, and 131 A channels along with the emissivities for the EIS high 
temperature lines (Fe XVI and Ca XlV-Ca XVII), and the two XRT channels considered in this 
paper. 

The AIA 94 A channel is of particular interest because it combines high spatial and temporal 
resolution with a relatively narrow temperature response. The EIS emission lines are formed at 
somewhat lower temperatures and, at least for these observations, have no temporal resolution. 
XRT observes these temperatures at high cadence, but has a very broad temperature response. 
Since the cross-calibration of the AIA instrument is just beginning at this point, we're interested 
only in some very basic questions, such as what is the morphology of the emission in this channel 
and how does it evolve with time? 

To provide some context for interpreting the AIA 94 A and 131 A channels we have used the 
simulations results for the high and low frequency heating cases shown in Figures [5] and [6] to 
compute the expected count rates in these channels. As one would expect, we find that the amount 
of high temperature emission observed in these channels is sensitive to the heating time scale. For 
the low-frequency heating case, which reaches temperatures close to 10 MK, observable emission 
(~ 50-100 DN s _1 ) is predicted for both channels. The high-frequency case shows some signal 
in the 94 A channel but relatively little in 131 A. These comparisons suggest that these channels 
will provide important information on the heating time scale. Quantitative comparisons, however, 
will require the AIA calibration to be more fully understood. Also, at these high temperatures the 
potential for departures from ionization equilibrium are much greater. 

It is clear from the movie of the 94 A channel, which is available in the electronic version of 
the manuscript, that there is significant emission from Fe XVIII in the core of the active region. 
This channel, however, also contains contributions from lines formed at lower temperatures. Com- 
parisons with the 171 A channel allow us to identify the high temperature, Fe XVIII loops. Such 
comparisons suggest that the bulk of the Fe XVIII emission lies at the inner edge of the moss on 
loops that connect directly across the neutral line. 

To address the question of temporal evolution we have taken the co-aligned and de-rotated data 
cubes used to make the movies and computed light curves for various points in the active region. 
Light curves for three points are shown in Figure [9j One of the points illustrates the evolution of 
a small flare. Here the light curves show the progression from the high temperature flare emission 
of Fe XX-Fe XXIII, to the hot Fe XVIII and Fe XVI, to the relatively cool Fe IX. This qualitative 
agreement between the observations and our expectations for a flare light curve suggests that the 
emission lines contributing to each channel are properly identified. Other small "microflaring" 
events show a similar progression, except that the cooling to the lowest temperatures is difficult to 
identify (middle panel). Finally, the majority of the pixels in the core of the active region do not 
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show coherent behavior. In these pixels the intensities measured in the 131 and 171 A channels are 
well correlated, suggesting an absence of plasma at flare temperatures. For these points the 94 and 
335 A channels do not show any obvious correlation. There is no general tendency for brightenings 
in the 94 A channel to be followed by brightenings in 335 A. This is potentially consistent with 
the high-frequency heating scenario, but will require more quantitative analysis to demonstrate 
conclusively. 



5. Discussion 

We have presented the detailed analysis of high temperature active region core loops using 
observations from Hinode and SDO. The primary result is the calculation of the emission measure, 
which shows a sharp peak at about 4 MK and a steep decline at both higher and lower temperatures. 
These observations provide a difficult challenge to the low-frequency heating scenarios that have 
traditionally been associated with nanoflare heating models of the corona. As we've shown with 
some very simple hydrodynamic simulations, such models predict significant emission at lower 
temperatures, but the observed emission measure below 1 MK is generally small in the core of an 
active region near the loop apexes. A high-frequency heating scenario, in contrast, agrees with the 
observed intensities of the high temperature emission lines, is consistent with the relatively steady 
emission observed in XRT, and with the absence of cooling loops in the moss in the cool AIA 
channels (e.g., 171 A). 

Our analysis provides compelling evidence that there is a significant population of loops in the 
core of a solar active region that are heated on very short time scales, much shorter than a typical 
cooling time for this combination of density (logn e ~ 9.7), temperature (logT e ~ 6.6), and loop 
length [L ~ 75 Mm) . We suggest that this combination of parameters is a useful set for theorists 
to consider in developing physical models of high temperature coronal loops. The challenge is 
to identify a physical mechanism that not only reproduces these conditions, but is also frequent 
enough to keep the loop at high temperatures for long periods of time. 

One limitation of the high-frequency heating scenario, however, is that it does not reproduce 
the emission that is observed at low temperatures. In the emission measure distributions shown 
in Figure for example, it is clear that the result from the high-frequency simulation is close to 
observed distribution at high temperatures (above logT = 6.5), but significantly under-predicts 
the distribution at low temperature s. One possib i lity is that there is a range o f heating time scales 



at work in the solar atmosphere. iBrooks et al.l (|2008l ) and iLee et al.l (|2010l ). for example, have 



suggested a relationship between magnetic topology and the variability of the heating. Transient 
brightenings appear to be associated with magnetic complexity, while magnetically simple regions 
show relatively constant emission. Motivated by the emission measure distributions shown in 
Figure H] we have constructed composite intensities with 90% high-frequency heating and 10% low- 
frequency heating. As is shown in Table this simple mixture provides a better match to the 
observations at both high and low temperatures. This doesn't provide a consistent description for 



-17- 



all of the lines, but does suggest how a complete description of the observed emission might be 
achieved. 

We stress that the hydrodynamic simulations presented here represent a very small set of pos- 
sible heating scenarios. Some OD impuls ive heating models have produced r elatively steep slopes in 



the emission measure distribution (e.g.. ICargilllll994l : iKlimchuk et al.ll2008l ). We have had trouble 



reproducing these results with hydrodynamic simulations and we are currently investigating this 
systematically. Other considerations, such as variations in loop length along the line of sight, differ- 
ences in initial conditions, or accounting for non-equilibrium ionization, may bring low-frequency, 
nanoflare heating models into better agreement with the observations. Much more extensive anal- 
ysis is required. One significant hurdle that all models must clear is reproducing the observed 
intensities both in the corona and at the footpoints. The next step in this analysis to couple 
hydrodynamic simulations to topology derived from magnetic field extrapolations and attempt to 
reproduce all of the emission observed in the active region. 

Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as domestic 
partner and NASA and STFC (UK) as international partners. It is operated by these agencies in 
co-operation with ESA and NSC (Norway). The authors would like to thank the referee for many 
valuable comments on the manuscript. 
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Table 1. Differential Emission Measure Model 1 
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a Calculated intensities are from the MCMC emission measure 
inversion. Intensities are in units of erg cm~ 2 s _1 sr -1 . 
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Fig. 6. — An example of low-frequency heating where the long time between heating events allows 
the loop to evolve over a wide range of temperatures. The top panels show the temperature and 
density averaged over the loop apex. The bottom panels show the evolution of the intensities in 
selected emission lines. The time averaged intensities, again taken over the final 8,000 s of the 
simulation are also indicated. The large intensities calculated for Fe XII 195.119 A and Si VII 
275.368 A are not consistent with the observations. 
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Table 2. Observed and Simulated Active Region Core Intensities' 
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a All intensities are in units of erg cm~ 2 s _1 sr _1 . The parameter R is the ratio of the observed to 
modeled intensities and r is time between heating events. Mixed refers to a combination of 90% high 
frequency and 10% low frequency intensities. 
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Fig. 7. — The emission measure distributions derived from the high- and low-frequency heating 
simulations. The distribution for both the entire loop (top panel) and loop apex (bottom panel) are 
shown. For comparison, our observed emission measure distribution is also shown in the bottom 
panel. The arrows indicate the differences between the observation and the low-frequency model 
at logT = 6.0 and 6.9. The power law indexes (EM ~ T b ) are indicated for several of the emission 
measure distributions. 
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Fig. 8. — Response curves for AIA, EIS, and XRT illustrating the temperature sensitivity of the 
various images and emission lines. The dotted vertical line indicates the peak temperature in the 
emission measure distribution. 
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Fig. 9. — Three representative light curves from the core of the active region. The top panels show 
the image data from the peak of the light curve. The bottom panels show the intensities as a function 
of time, the intensities have been smoothed with 60s running averages. The first column shows 
the evolution of a small flare, with a clear progression from high temperatures (Fe XX-Fe XXIII) to 
low temperatures (Fe IX). The middle column shows a similar behavior for a microflare. The final 
column shows the time history of a relatively quiescent pixel in the active region core. 



